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Abstract 



Bushes of normal modes represent exact mathematical objects describing specific 
dynamical regimes in nonlinear mechanical systems with point or space symmetry. 
^ ' In the present paper, we outline the bush theory and illustrate it with some bushes 

of small dimensions in octahedral mechanical structures. 



1 Introduction 



1. 



A new concept, "bushes of normal modes," was introduced for nonlinear mechanical sys- 
tems with discrete symmetry in ^ . A given bush represents a certain superposition of 
I i: the modes associated with different irreducible representations (irreps) of the symmetry 
! group G of the mechanical system in equilibrium. The coefficients of this superposition 
Q I are time-dependent functions for which the exact ordinary differential equations can be 
^ ' obtained. In this sense, the bush can be considered as a dynamical object whose dimen- 
' sionality is frequently less than that of the original mechanical system. The following 
^ . propositions were justified in previous papers [|l], ^ ^: 

A certain subgroup Gd of the symmetry group G corresponds to a given bush, and 
this bush can be excited by imposing the appropriate initial conditions with the 
above symmetry group G D a G. 

Each mode belonging to the bush possesses its own symmetry group which is greater 
than or equal to the group Gr, oi the whole bush. 

In spite of evolving mode amplitudes, the complete collection of modes in the given 
bush is preserved in time and, in this sense, the bush can be considered as a geo- 
metrical object. 

The energy of the initial excitation is trapped in the bush, i.e. it cannot spread to 
the modes which do not belong to the bush, because of the symmetry restrictions. 

As an indivisible nonlinear object, the bush exists because of force interactions 
between the modes contained in it. 



6. 



Taking into account the concrete type of interactions between particles of the con- 
sidered mechanical system can only reduce the dimension of the given bush. 



7. The extension of the bush can be reahzed as a result of the loss of its stability which 
is accompanied by spontaneous breaking of the bush symmetry (dynamical analog 
of phase transition). 

The special group-theoretical methods for finding bushes of modes are discussed in [|, 
^. The computer implementation of these methods (see also 0) allowed us to find 
bushes of modes for wide classes of mechanical systems with discrete symmetry. In par- 
ticular, "irreducible" bushes of vibrational modes and symmetry determined similar non- 
linear normal modes for all N-particle mechanical systems with the symmetry of any of 
the 230 space groups were found in IQ. The bushes of vibrational modes of small dimen- 
sions were found and classified into universality classes for all mechanical systems with 
point groups of crystallographic symmetry in 0. 



2 Outline of the bush theory 

We consider classical Hamiltonian systems of mass points moving near the single equi- 
librium state which can be characterized by a certain point or space symmetry group G. 
Let the 3 x A^-dimensional vector, 

X(t) = (Xi(t),X2(t),...,X^(t)), (1) 

describe the displacements Xj(t) of all particles of our mechanical system from their equi- 
librium positions. (Here we denote by the three-dimensional vector Xj the displacement 
of the i-th particle along the X, Y and Z axes). 

The vector X(if:) can be written as a superposition of all basis vectors (p^P of the 
irreducible representations Tj of the above mentioned symmetry group G: 

x{t)=j:f^^it}^?- (2) 

The coefficients ^^i\t) of this superposition depend on time t, while the 3 x A^-dimensional 
time-independent vectors c^) are determined the specific patterns of displacements of all 
particles of our mechanical system. 

Note that individual components of the basis vectors ^pf^ are often called "symmetry- 
adapted coordinates" . In particular, they can be normal coordinates. Hereafter, the term, 
"mode," means an arbitrary superposition of basis vectors corresponding to a given irrep 
Tj. As a result of this definition, every term jji^p{t)ip^p in the right hand side of Eq.(|]) is 
also a mode of the irrep Tj. Sometimes, for brevity, we will refer to ^^i\t) as a mode, but a 
reader must imagine that this time-dependent coefficient is multiplied by the appropriate 

3 X A^- dimensional vector to give the mode in the exact sense. We can also speak 
about vibrational modes because only such type of symmetry-adapted (normal) modes 
are considered in the present paper. 

Every dynamical regime of the considered mechanical system can be described by 
the appropriate time- dependent vector X(t) which determines a definite instantaneous 
configuration of the system. On the other hand, each instantaneous configuration pos- 
sesses a certain symmetry group Go (in particular, this group may be trivial: Gr, = 1) 
which is a subgroup of the symmetry group G of the system in equilibrium {Gd C G). 
Moreover, we can also ascribe a certain symmetry group to each basis vector cp'f'^ and to 
each mode corresponding to a given irrep Tj (remember that a mode is a superposition 



of such vectors!), because the definite instantaneous configurations correspond to them. 
The group Gr, contains all symmetry elements of group G whose action does not change 
this configuration. 

Let us introduce, as it is usual in group theory, the operators g associated with elements 
g of group G {g E G) which act on 3 x A^-dimensional vectors X(t).Q All elements g E G 
for which 

gX{t) = X(t) (3) 

form a certain subgroup Gd <^ G, and a complete set of the above operators g (Wg E Go) 
represents the group Gd- 

It can be shown that the symmetry group G/j is preserved in time in the sense that its 
elements cannot disappear during time evolution. Actually, this property is a consequence 
of the principle of determinism in classical mechanics.0 Thus, the equation gX{t) = X(t) 
{g G Gd), or, formally, 

GoMt) = X(t) (4) 

is valid for every time t. As a consequence, we can classify the different dynamical regimes 
in our nonlinear dynamical system, described by the vectors X(t) from Eqs.(|I|,0), with 
the aid of symmetry groups corresponding to them. 

Using Eq.@, we can obtain the similar mvariance conditions for each individual irrep 
Vj of the group G (see details in 0): 

(rj i G'd)Mj = Mi- (5) 

Here {Vj [ Gd) is a restriction of the irrep Vj to the subgroup Gd of the group G, and 
, /i^^,)) is an invariant vector oi Vj {uj is the dimension of this irrep). 

To find all modes contributing to a given dynamical regime with symmetry group 
Gd, i.e., for the vector X(t) from Eq.(0), we must solve linear algebraic equations 
for each irrep Vj of the group G. As a result of this procedure, the invariant vector fij 
for some irreps Vj can turn out to be equal to zero. Such irreps do not contribute to 
the considered dynamical regime. On the other hand, some nonzero invariant vectors fij 
for multidimensional irreps may be of a very specific form because of definite relations 
between their components (for example, certain components can be equal to each other, 
or differ only by sign). 

Actually, Eq.(|^) can be considered as a source of certain selection rules for spreading 
excitation from the root mode to a number of other (secondary) modes. Indeed, if a 
certain mode with the symmetry group G^ is excited at the initial instant (we call it the 
"root" mode), this group determines the symmetry of the whole bush. The condition that 
the appropriate dynamical regime X(t) must be invariant under the action of the above 
group Gd leads to Eq.(^ and then to Eq.(|^). If the vector /x^- for a given irrep Vj proves 
to be a zero vector, then there are no modes belonging to this irrep which contribute to 
X(t), i.e., the initial excitation cannot spread from the root mode to the secondary modes 
associated with the irrep Vj. 

Note that basis vectors associated with a given irrep Vj in Eq.(|^) turn out to be equal 
to zero when this irrep is not contained in the decomposition of the full vibrational irrep 
r into its irreducible parts Vj. This is a source of the additional selection rules which 
reduce the number of possible vibrational modes in the considered bush. Trying every 

"'^The symmetry elements g € G act on the vectors of three-dimensional Euclidean space. 
^The phenomenon of spontaneous breaking of symmetry of a given dynamical regime will be considered 
in the next section. 



irrep Tj in Eq.(|^) and analyzing the above mentioned decomposition of the vibrational 
representation F, we obtain the whole bush of modes with the symmetry group Go in the 
explicit form. 

Let us return to Eq.(|D. We speak about geometrical aspects of the bush theory when 
concentrate our attention on basis vectors iPi'\ and we speak about dynamical aspects of 
this theory when we focus on time-dependent coefficients fi\''\t), which will be also called 
"modes" . 

If interactions between the particles of our mechanical system are known, exact dy- 
namical equations describing the time evolution of a given bush can be written. 

Two types of interactions between modes in nonlinear Hamiltonian system are dis- 
cussed in 1^, namely, force interactions and parametric interactions. We can illustrate the 
difference between these types of modal interactions using a simple example. 

Let us consider two different linear oscillators whose coupling is described by only one 
anharmonic term, U = —7/2^/12, in the potential energy. Dynamical equations for this 
system can be written as follows: 

/ii+cj^/ii = 2jfiiti2, (6) 

il2+i^lfi2 = llA- (7) 

Here 7 is an arbitrary constant characterizing the strength of the interaction of the os- 
cillators. We can suppose that Eqs.(|],^ describe the dynamics of two modes /ii(t) and 
/i2(t) in a certain mechanical system. 

An essential disparity between modes and ^2{t) can be seen from the above 

equations. Indeed, if we excite the mode ^i{t) at the initial instant (/ii(to) 7^ 0), the 
mode ^2if) cannot be equal to zero (even if it was zero ait = t^l) because a nonzero force 
~ ^''^1 appears in the right hand side of Eq.(^ since ^ 0. In other words, 

the dynamical regime ^i{t) ^ 0,/i2(t) = cannot exist because of the contradiction with 
Eq.(^. Unlike this, the dynamical regime fi2(t) ^ 0,/ii(t) = can exist because such 
a condition does not contradict equations (|^) and (0). We can say that now there is no 
force in the right hand side of Eq.(0) because ^ = as a consequence of the identity 
^i(t) = 0. 

In the last case, the dynamical regime of the system represents a harmonic 

oscillation only of the second variable: 

H2it) = Acos{uj2t + 5), (8) 

where A and 6 are two arbitrary constants. 

Thus, there is force action from the mode Hi(t) on the mode /U2(t), but not vice versa. 
We proved in that such a situation can be realized only in the case where the symmetry 
group of the mode fii{t) is less than or equal to that of the mode /i2(^)-0 

Nevertheless, the mode fi2{t) can excite the mode fii{t) under certain circumstances. 
Indeed, substituting the solution of Eq.(|^ into Eq.(^, we obtain 

ili{t) + \ul - 2yl7cos(cj2t + S)] /Ui = 0. (9) 



■^Remember, that speaking about a mode we must take into account that our time-dependent coefficient 
is multipHed by the appropriate basis vector of a certain irrep of the group G, and namely this vector 
determines the symmetry group of the considered mode. 



By means of simple algebraic transformations this equation can be converted to the Math- 
ieu equation in its standard form: 

+[a-2q cos(2r)] z = 0, (10) 

where z = z{t). But in the (a — q) plane of the Mathieu equation (|l^) there exist domains 
of stable and unstable movement. If pertinent parameters [ui,U2,'y) of our dynamical 
system have values such that corresponding parameters a and q of Eq. (|lOD get into an 
unstable domain, then the nonzero function z{t) and, therefore, the mode fii{t) appears. 
In other words, the initial dynamical regime ^i{t) = 0, //2(i) ^ loses its stability, and a 
new dynamical regime ^ 0,/i2(t) ^ arise spontaneously for definite values of the 

parameters of Eqs.(^|^. Since this phenomenon is similar to the well-known parametric 
resonance, we can speak, in such a case, about parametric action from the mode /X2(^) on 
the mode 

The characteristic property of the parametric interaction is that the appropriate force 
(— ^ = 27/ii/i2, in our case) vanishes when the mode (//i, in our case), on which this 
force acts, becomes zero. The following important result was proved in the mode 
of lower symmetry acts on the mode of higher symmetry by force interaction, while the 
mode of higher symmetry can act on the mode of lower symmetry only parametrically. 
Consequently, if the parametric excitation of a certain mode does take place, this phe- 
nomenon must be by necessity be accompanied by spontaneous breaking of symmetry of 
the mechanical system vibrational state. 

Thus, the initially excited dynamical regime can lose its stability because of para- 
metric interactions with some zero modes and, as a result, can transform spontaneously 
into another dynamical regime, described by a greater number of dynamical variables, 
with appropriate lowering of symmetry. Obviously, we may treat such phenomenon as a 
dynamical analog of a phase transition. 

3 Examples of bushes of vibrational modes 

We consider a mechanical system of six mass points (particles) whose interactions are 
described by a pair isotropic potential u{r) where r is the distance between two particles. 
We suppose that in the equilibrium state these particles form a regular octahedron with 
edge Oo which can be imagined in the following way. Let us introduce a Cartesian coor- 
dinate system. Four particles of the above octahedron lie in the XY plane and form a 
square with edge a^. Two other particles lie on Z axis and we will speak about the "top 
particle" and the "bottom particle" with respect to the direction of the Z axis. Obviously, 
the distance between each of these two particles and any of the four particles in the XY 
plane is equal to oq. 

The point symmetry group of the octahedral equilibrium configuration is Oh and 
all bushes of vibrational modes in the considered system are described by the certain 
subgroups Go oi this parent group. In the present paper, we consider only three bushes^ 
Bl[Oh\, B2[D/ih\ and i?4[C4^] (the complete list includes 18 different by symmetry bushes 
of vibrational modes whose dimensions vary from 1 to 12). The symmetry groups of the 
above bushes are connected with each other by the following group- subgroup relation: 
6*4^ C -D4/1 C Oh- 

^We write the symmetry group Gd of the bush in square brackets next to its symbol. 



The geometrical forms of our mechanical system in the vibrational state, corresponding 
to these bushes, can be revealed from the appropriate symmetry groups Ge>. 

The one-dimensional bush i?l[0/j] consists of only one ("breathing") mode. This 
nonlinear dynamical regime X(t) = {t)(pi^^ describes evolution of a regular octahedron 
whose edges a — a{t) periodically change in time. 

The two-dimensional bush B2[D4h] describes a dynamical regime with two degrees 

of freedom: X(t) = A*i^'*(^)¥'i^'* + A*i^'*(^)¥'i^''- The symmetry group Gd = D^h of this 
bush contains the 4-fold axis coinciding with the Z coordinate axis and the mirror plane 
coinciding with the XY plane. This symmetry group restricts essentially the form of 
the polyhedron describing our mechanical system in the vibrational state. Indeed, the 
presence of the 4-fold axis demands that the quadrangle in the XY plane be a square. 
Because of the same reason, the four edges connecting the particles in the XY plane 
(vertices of the above square) with the top particle lying on Z axis must be of the same 
length which we denote by b{t). 

Similarly, let the length of the edges connecting the bottom particle on the Z axis 
with any of the 4 particles in the XY plane be denoted by c{t) . In our present case of the 
bush B2[D4h], b{t) = c{t) for any time t because of the presence of the horizontal mirror 
plane in the group Gu = D^^. But for the three-dimensional bush i?4[C4t,], described by 
X(t) = iJ'i^\t)(p^i^ + lJ'f\t)(pf^ + iJ,^^^\t)(f^^\ this mirror plane is absent and, therefore, 

b{t)y^cit). 

Let us also introduce two heights, hi{t) and h2{t), corresponding to the perpendiculars 
dropped, respectively, from the top and bottom vertices of our polyhedron in the XY 
plane. Now we can write the dynamical equations of the above bushes in terms of pure 
geometrical variables a{t) , b{t) , c{t) , hi{t) and h2{t). 

We choose a{t) and h{t) = hi{t) = h2{t) as dynamical variables for describing the two- 
dimensional bush B2[D4^h\ and a(t), hi{t) and h2{t) as dynamical variables for describing 
the three-dimensional bush i?4[C4^] . Using these variables we can write down the potential 
energy for our bushes of vibrational modes as follows: 

Bl[Oh] : VBi(a) = 12M(a) + 3ii(v^a), 

B2[D4h\ ■■ yB2(a,/i) = 4ii(a) + 2M(V2a)+8M(y^^M^) +^x(2/i), 
BA[Ci^] : VB4(a,/ii,/i2) =4M(a) + 2M(^/2a)+4^6)+4ii(c) + M(/ii + /i2), 

where b = yfT^f^^T^i^, c = + (|/i2 - i/^i)'- 

Then with the aid of the Lagrange method, we can obtain the following dynamical 
equations for the above bushes of vibrational modes: 



Bl[Oh] : 






d — 


-Au'{a) - 


\/2M'(\/2a); 


B2[D4f,] : 






d — 


-2u'(a) - 
-4«'(6)f- 


V2u'{^/2a) -2u'(b)^, 


h = 


-u'{2h)] 


54[C4„] : 




d — 


-2u'{a) - 


V2u'{V2a) ^u'{b)l - 




-u'ib)'-^ 


^-U'{hi + h2), 


h2^ 


m'(c)5'^2- 


-U'{hi + h2). 



Thus, we obtain the dynamical equations of our bushes of vibrational modes in terms 
of variables with explicit geometrical sense. Each bush describes a certain nonlinear 



dynamical regime corresponding to such a vibrational state of the considered mechanical 
system, that at any fixed time the configuration of this system is represented by a definite 
polyhedron with symmetry group Gr, of a given bush. 

We can write dynamical equations for the above bushes in terms of vibrational modes 
as well. In spite of the more complicated form, these equations turn out to be more 
useful for the bush theory, since they allow us to decompose the appropriate nonlinear 
dynamical regimes into modes of different importance for the case of small oscillations - 
root modes and secondary modes of different orders 0. As an example, we write below 
the dynamical equations for the bush S4[C4i;] in terms of its three modes, n^^^\t) = 7(t) 
(root mode), fJ'^i\t) = fJ'i'\t) = iy{t) (secondary modes): 

il = -|(4v^M'(a) + 4u'{V2a) + 2u'{hi + + ^(2v^a + 5h - /i2) + 

^{2V2a + 5h2-h)), 
i) = -|(2\/2u'(a) + 2u'{V2a) - 2u'{hi + /is) + ^(v^a - 5hi + h2) + 

Here 

a = V2{ro + + u);b = ^{rp + fi + uf + {rp + fi - 2iy - 87)^; 
c= y^(ro + // + z/)2 + (ro + /i-2// + 37)2; 
hi = ro + fi-2u-2-f;h2 = ro + fi-2u + 27; Vp = ^. 

In the previous section, the problem of existence of bushes of vibrational modes was 
studied by group-theoretical methods only. In contrast to this, examination of the stability 
of the bushes depends essentially on the concrete type of interactions in the considered 
system, and we suppose that they can be described by a Lennard- Jones potential. Then 
the bush stability will be analyzed with the aid of numerical methods. 

Let us excite a given bush B[Gd] with the aid of the following initial condition. Note 
that the root mode of each of the above considered bushes is determined by a single time- 
dependent coefficient ^p'' (t) whose initial value at t = tp we will denote by the symbol ^p 
{/j,p = fi[-'\tp)). At the initial instant t = tp, we fix the coordinates of all particles of the 
mechanical system in such a way that their displacements correspond to the appropriate 
root mode with amplitude fip, while their velocities are equal to zero. Namely this choice 
of initial conditions determines the way of excitation of a given bush. 

Using these initial conditions we solve numerically the exact dynamical equations of 
the considered mechanical system with 18 degrees of freedom, and analyze the set of 
nonzero modes ^i'\t) in the decomposition (^ of the vector X(t) obtained as a result 
of this solving. Then we gradually increase the value fip and repeat the procedure just 
described until the number of nonzero modes ^[''\t), at some value /xq = R, becomes 
larger than that of the bush B[Gj:)]. We will refer to R as the threshold of stability of the 
given bush. Obviously, in such a way we obtain the upper boundary of the first stability 
region of B[Gd] in one-dimensional space of all possible values fip (0 < /io < R)-^ 

^Note that we do not study the other possible regions of stabihty of the given bush B[Gd] in the 
present paper. 



For fiQ > R a new bush B[G£)], including the old bush B[G£)], appears and its symme- 
try Ge> is lower than that of this old bush {Gu C Gn), because all modes with symmetry 
higher than or equal to Gd are already contained in B[Gd]. 

Thus, the loss of stability of a given bush is accompanied by the spontaneous breaking 
of symmetry of the initially excited dynamical regime, described by this bush. We already 
discussed this phenomenon in Sec.|] and concluded that its cause is analogous to that of 
this parametric resonance. 

We can also say this in other words. A given bush B[Gd] represents a certain dynam- 
ical regime in the considered mechanical system. Its modes interact with other modes 
which do not belong to B[Gd], but these interactions must be of parametric (not force!) 
type only. For the appropriate initial conditions we can get into a region of unstable 
movement. As a result, some new modes are excited which were forbidden by principle of 
determinism of classical mechanics. Then we can speak about the loss of stability of the 
original bush B[Gd] and its transformation into a larger bush B[Gd] with {Gd C Gd)- 

The Lennard- Jones potential can be written in the form, 

= —2--,^ (12) 

with A = B = 1 (because of the appropriate scaling transformation of space and time 
variables). Our numerical experiments give the following values of the threshold R for 
the above considered bushes: 



R[B1] = 0.001, R[B2] = 0.009, i?[54] = 0.003. (13) 

Note that these values of R are finite but very small in comparison with the edge 
ao = 1.117 of our octahedral structure in equilibrium. Displacements of particles corre- 
sponding to the thresholds (|T^ can be obtained by dividing these three values of R by 
the numbers ^/6, a/T2, vT2, respectively. We want to note in this connection, that all 
octahedral molecules, known to us at the present time, possess an atom in the center of 
the octahedron. This fact suggests that the stability of such structures can be greater 
than that of the mechanical system considered up to this point. Taking into account this 
hint, we examined bush stability for the mechanical structure with the particle in the 
center of the octahedron, supposing that this additional (seventh) particle is described by 
the Lennard- Jones potential different from that for six peripheral particles. Namely, we 
assume that A = 1, but B > 1 in Eq. (p!2D . Such an assumption provides us a possibility 
of making the attractive part of the potential of the centered particle greater than that 
of peripheral atoms in spite of the same repulsive part. 

For such a centered structure the threshold values R can be essentially greater than 
those for the structure without the particle in the center of the octahedron. Indeed, we 
obtained that R[B1] = 1.010, R[B2] = 0.011, R[BA] = 0.118 for B = 5.5. Note that the 
dependence of the threshold R on the value B is essentially different for different bushes 
and can be nonmonotonic. 
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